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In this paper we present an algorithm to reduce the area of a surface spanned by a finite number 
of boundary curves by initiating a variational improvement in the surface. The ansatz we suggest 
consists of original surface plus a variational parameter t multiplying the numerator Ho of mean 
curvature function defined over the surface. We apply this technique to a hemiellipsoid and bilinear 
interpolation spanned by four bounding straight lines. (The four boundary lines of the bilinear inter- 
polation can model the initial and final configurations of re-arranging strings.) As a demonstration 
of the effectiveness of the technique, the area of the hemiellipsoid is reduced for the same boundary 
by as much as 23 percent of the original area. For bilinear interpolation the decrease remains less 
than 0.8 percent of the original area, which may suggest that it is already a near minimal surface. 



I. INTRODUCTION 

Variational methods are one of the active research areas of the optimization theory [l| . A variational method tries 
to find the best values of the parameters in a trial function that optimize, subject to some algebraic, integral or 
differential constraints, a quantity dependant on the ansatz. A simple example of such a problem may be to find 
the curve of shortest length connecting two points. The solution is a straight line between the points in case of no 
constraints, otherwise possibly many solutions may exist depending on the nature of constraints. Such solutions are 
called geodesies. One of the related problems is finding the path of stationary optical length connecting two points, as 
the Fermat's principle says that the rays of light traverse such a path. Another related problem is a Plateau problem 
0] which is finding the surface with minimal area enclosed by a given curve. This problem is named after the blind 
Belgian physicist Joseph Plateau, who demonstrated in 1849 that a minimal surface can be obtained by immersing 
a wire frame, representing the boundaries, into soapy water. The Plateau problem attracted mathematicians like 
Schwarz Q (who discovered D (diamond), P (primitive), H (hexagonal) , T (tetragonal) and CLP (crossed layers of 
parallels) triply periodic surfaces), Riemann and Weierstrass 0. Although mathematical solutions for specific 
boundaries had been obtained for years, but it was not until 1931 that the American mathematician Jesse Douglas Q 
and the Hungarian Tibor Rado [|| independently proved the existence of a minimal solution for a given simple closed 
curve. Their methods were quite different. Douglas Q minimized a quantity now named as Douglas Integral and his 
results held for arbitrary simple closed curve, while Rado [f| minimized the energy. The work of Rado was built on 
the previous work of R. Gamier Q and held only for rectifiable simple closed curves. Many results were obtained in 
subsequent years, including revolutionary achievements of L. Tonelli i, R. Courant 01 [lO] , C- B. Morrey [ll[ 
E. M. McShane [3, M. Shiffman [3, M. Morse [H], T. Tompkins [lj|, Osserman [lj|, Gulliver [H and Karcher 
[3 and others. 

An emerging use of minimal surfaces in physics is that in string theories. A classical particle travels a geodesic 
with least distance whereas a classical string is an entity which traverses a minimal area. Amongst the string theories 
used in physics, two are worth mentioning. One is the theory of quantum chromodynamics (QCD) strings that model 
the gluonic field confining a quark and an antiquark within a meson. (The gluonic field connecting three quarks, 
within a proton or neutron, is modeled through Y-shaped strings. For a system composed of more than three quarks, 
minimization of the total length of a string network with only Y-shaped junctions may be a non-trivial Steiner-Tree 
Problem [T9|). In the other string theory (or theories) string vibrations are supposed to generate different elementary 
particles of the present high energy physics. Quite often string theories need a surface spanning the boundary composed 
of curves either connecting particles or describing the time evolution of particles. An important case can be a fixed 
boundary composed of four external curves. A common application of this boundary can be the time evolution of a 
string parameterized by a or (3 [2(| variable; the time evolution itself is parameterized by the symbol r, the proper 
time of relativity. In this case two bounding curves parameterized by the respective a or j3 represent the initial and 
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final configurations of a string, and the other two curves (parameterized by the respective r variables) describe the 
time evolution of the two ends of a string. 

String theories take action to be proportional to area. Combining this with the classical mechanics demand of 
the least action, minimal surfaces spanning the corresponding fixed boundaries get their importance. For example, 
see eq. 13 of ref. [20[ for the Nambu-Goto ansatz for the minimal surface area and compare it with eqs. (fT2"]) and 
(fl~3|) below, along with ref. [2l| for Nambu-Goto strings. Also relevant is the use in ref. |22j of Wilson minimal 
area law (MAL) to derive the quark antiquark potential in a certain approximation. A surface spanned by such a 
boundary is in space-time of relativity. An ordinary 3-dimensional spatial surface can span a boundary composed 
of two 3-dimcnsional curves connecting four particles and two other curves connecting the same four particles in a 
re-arranged (or exchanged) clustering; see for example Fig. 2 of ref. [23[ and Fig. 5 of ref. [24| . An explicit expression 
of such a spanning surface can be found in eqs. 3, 4 of ref. J2(| and eq. 22 of ref. (25|. This is a bilinear interpolation 
in ordinary 3-dimensional space and is similar to the linear interpolations in above mentioned eq. 13 of ref. [20| . 
eq. 4.7 of ref. [28[ and eq. 3.4 of ref. [22j|. R- ei - 113 clarifies that such a surface is used as a replacement to the 
exact minimal surfaces for the corresponding boundaries; see section [TT] below for a minimal surface in the differential 
geometry. Even non-minimal surfaces have some usage in the mathematical modeling of quantum strings because 1) 
in contrast to classical strings, quantum strings can have any action and hence area as described by the path integral 
version of the quantum mechanics (see eq. 1 of [3~l| ) and 2) any surface spanning a boundary composed of quark lines 
(or quark connecting lines) corresponds to a physically allowed (gauge invariant) configuration of the gluonic field 
between these quarks; compare the non-minimal surface of Fig. 10.5 of ref. (29| with the minimal surface for the same 
boundary in Fig. 10.1 of the same ref. (29|. But it cannot be denied that minimal surfaces are the most important of 
the spanning surfaces even in quantum theories. For example, the relation in eq. 1.14 of ref. [3(| between an area and 
an important quantity (termed Wilson loop) related to the potential between a quark and antiquark connected by a 
QCD (gluonic) string is valid only if the area is of the minimal surface. (Thouth above mentioned eq. 1 of ref. [3l| 
relates the Wilson loop to a "sum over all surfaces of the topology of rectangle bounded by the loop" implying that 
each spanning surfaces has some contribution in the Wilson loop, the minimal surface must contribute most.) Thus it 
is worth pointing out that the non-minimal linearly or bilinearly interpolating surfaces can replace minimal surfaces, 
can be effectively used as minimal surfaces or share some features in common with minimal surfaces; text just before 
eq. (1.15) of the above mentioned ref. (30j relates them, up to non-relativistic l/(mass square) order, to the minimal 
surfaces. The purpose of the present paper is explore further this "effective usability" or "sharing common features 
with minimal surfaces" of linearly or bilinearly interpolating surfaces. Before starting a description of our work, we 
want to 1) state the common feature we have chosen. This is the fractional reduction possible in the area for a fixed 
boundary; for an exact minimal surface this quantity is zero (at least for a small neighbourhood). For reducing area 
we use the variational area reduction, outlined in sect. IIV1 to our specific bilinear interpolation described in sec. IIIII 
Moreover, we 2) point out that the bilinear interpolations used in string-theories-related works of physics are also used 
in the emerging discipline of the computer aided geometric design (CAGD) and hence the usefulness of the present 
paper extends to CAGD along with physics and the differential geometry; as much as bilinear interpolations are near 
or related to minimal surfaces their study sheds some light on the above mentioned Plateau problem of the differential 
geometry itself. 

Computer aided geometric design (CAGD) [Hj], [33| arose when mathematical descriptions of shapes facilitated the 
use of computers to process data and analyze related information. In the 1960s, it became possible to use computer 
control for basic and detailed design enabling utilization of a mathematical model stored in a computer instead of 
the conventional design based on drawings. The term geometric modeling is used to characterize the methods used 
in describing the geometry of an object. Over the years, various schemes were developed with a view to achieve this 
abstraction. S. A. Coons (32|, [33| introduced the Coons patch in 1964. The Coons patch approach is based on the 
premise that a patch can be described in terms of four distinct boundary curves. Thus a Coons patch can be a worth 
analyzing surface spanning a fixed set of boundary curves, at least when the number of curves is four. For us, Coons 
patch (see eq. (|9]) below) is relevant because the above mentioned bilinear interpolations (see also eq. (fTOf below) we 
basically study in this paper are a special case of Coons patch [32[ . Coons patch analysis is an active area of research 
and has seen enormous development during recent years. But most, if not all, of the work on it has been limited 
to its geometric descriptions and visualization and to interactive mathematical experiments with it; it has not been 
analyzed from the view of differential geometry and that is also what we aim to do in this paper though we actually 
study only its special case of a bilinear interpolation. 

The paper is organized as follows. In sections [TT] and IIIII we present basic definitions and constructions related to 
surfaces spanned by fixed boundary curves. In the next section IIVI we present an algorithm to reduce the area of a 
surface spanned by a finite number of boundary curves by introducing a variational improvement in a surface. Then 
in section [V] we apply this technique to reduce the area of a non-minimal surface spanning a boundary for which 
the minimal surface is known - namely hemiellipsoid eq. (|28|) , to make sure the efficiency of the algorithm given 
by eq. (f2"Tj) and above mentioned bilinear interpolation spanned by four bounding lines for which the corresponding 



3 



minimal surface is not known. Based on this comparison, we comment on the possible status of bilinear interpolation 
as an approximate minimal surface. The last section IVTI presents results, final remarks and mentions possible future 
developments. 



II. DIFFERENTIAL GEOMETRY OF MINIMAL SURFACES 



In the optimization problem we aim for here, we eventually try to find a surface of a known boundary that has a 
least value of area. Area is evaluated by the area functional: 

A(x) = J J \x u (u,v) x x v (u,v)\ dudv, (1) 

where D C R 2 is a domain over which the surface x(u, v) is defined as a map, with the boundary condition x(dD) = T 
for < u < 1 and < v < 1, x u (u, v) and x v (u, v) being partial derivatives of x(u, v) with respect to u and v . It is 
known [35j~that the first variation of A(x) vanishes everywhere if and only if the mean curvature H of x(u, v) is zero 
everywhere in it. Thus a surface of least area is also a surface of least (zero) rms mean curvature spanning the given 
boundary. This means we can aim for the same surface using the condition of the least means square mean curvature 
in place of the condition of the least area. This is helpful as, unlike area, the ms mean curvature has not a square 
root in its integrand. For a locally parameterized surface x = x(u, v), the mean curvature H may be given by 

Ge-2Ff + Eg 

H ~ EG - F 2 ' (2) 

where 

E=(x„,x„), F=(x„,x„), G = (x v ,x v ), (3) 
are the first fundamental coefficients and 

e = (N, x uu ) , / = (N, x uv ) , g = (N, x„) , (4) 
are the second fundamental coefficients with 

being the unit normal to the surface x(u, v). For a minimal surface [35j . [36j j the mean curvature ([2]) is identically zero. 
For minimization we use only the numerator part of mean curvature H given by ((2|), as done in rcf. [37| following 
ref. [H who writes that "for a locally parameterized surface, the mean curvature vanishes when the numerator part 
of the mean curvature is equal to zero" . We call this numerator part Hq as the rms curvature of the initial surface 
x (w, v) to be used in the ansatz eq. (|2Tj) to get first order variationally improved surface Xi(w, v) of lesser area. This 
process could be continued as an iterative process until a minimal surface is achieved. But due to complexity of the 
calculations required for obtaining the second order improvement X2(w, v), we have been able to calculate the first 
order Xi (u, v) only. The numerator part Hq is denoted by 

H = e G - 2F f +g E , (6) 

where Eq, Fq, Go, eo, /o and go denote the fundamental magnitudes given by eqs. ([3|) and (j4|), with Nq(u,v) being 
the unit normal given by eq. §5§ to the initial surface x (it,i>). We call the root mean square (rms) of this H , for 
< u < 1 and < v < 1, as fj,Q. That is, 

1/2 

/ / / TT 9 , , ^ 



o Jo 



Hi dudvj . (7) 

In the notation of eqs ([2]) to (J4j eq. (pj becomes, for Xo(it, w), 

A = ( f Je q G -F 2 dudv. (8) 
Jo Jo v 
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III. BILINEAR STARTING SURFACE SPANNED BY FIXED BOUNDARY CURVES 



For a minimal (or, more precisely, a stationary) surface, we have to solve the differential equation obtained by setting 
the mean curvature H given by eq. @ equal to zero for each value of the two parameters, say, u and v parameterizing 
a surface spanning the fixed boundary. In this section, our purpose is to describe a starting surface bounded by the 
skew quadrilateral which is composed of four arbitrary straight lines connecting four corners x(0,0),x(0, l),x(l,0) 
and x(l,l); in the next section we report the variational improvement to this start aimed towards minimizing the 
surface evolving from the starting surface of this section. Ref. p7j also includes a preliminary effort to variationally 
improve the surface bounded by four straight lines towards being a minimal surface. The algorithm used for this 
variational improvement applies to a wider class of surfaces. Accordingly, now we point out a class of surfaces, namely 
Coons patch, that includes surfaces bounded by four straight lines: Let ci(m), c 2 (u), di(v) and d2(i>) be four given 
arbitrary curves defined over the parameters u,v S [0, 1] . For ci(u) = x(u,0), c 2 {u) = x(u, 1), di(v) = x(0, u) and 
d- 2 {v) = x(l,u), blending functions fx(u), f 2 (u), gi(v) and g 2 (v) satisfying the conditions that fx(u) + f 2 (u) = 1, 
gi(v) + g 2 (v) = 1 for non-barycentric combination of points and /i(0) = gi(0) = l,/i(l) = gi(l) = in order to 
actually interpolate x(0, 0), x(0, 1), x(l, 0) and x(l, 1), the following equation defines Coons patch: 



l{u,v) = [ fi(u) f 2 (u) 

- [/i(u) f 2 (u) 



x(0,v) 
x(l,w) 

x(0,0) x(0,l) 

x(l,0) x(l,l) 



+ [ x(u,0) x(w, 1) 



9i(v) 
92(v) 



9i{v) 
92{v) 



(9) 



As a special case of the above, we consider a Coons patch for which all the three terms in eq. © are equal, so that 
this equation reduces to the following form: 



x(u,v) = [ fi(u) f 2 (u) 



x(0,0) x(0,l) 
x(l,0) x(l,l) 



9i(v) 
92(v) 



(10) 



The boundary spanned by lines connecting the points x(0, 0), x(0,l), x(l,0) and x(l,l) with linear blending 
functions 



/l = 1 - u, f 2 =u, gi = l- v, g 2 = v, 



(11) 



in eq. (|10p can represent a time evolution of a string or, alternatively, a re-arrangement of a one set of two strings 
to the only possible other re-arranged set (of two strings) connecting the same two particles and two antiparticles. 
(It is to be noted that a string connects only a particle with antiparticle. This constraint allows only two string 
arrangements for a system composed of two particles and two antiparticles.) Above x(u, v) eq. (fTU|) spanning this 
boundary is a surface that is needed in many models of string re-arrangements from one of these configurations to the 
other with the particle positions x(0, 0) = 1 and x(l, 1) = 2 and anti- particle positions x(l,0) = 3 and x(0, 1) = 4. 
In this paper we reduce the the area of a quadrilateral, using above linear blending functions and particle positions. 
This gives 



Xn(u,v) = (1 - v)r 1 



vr 2 z, 



and 



x v (u, v) = (1 - u) r l3 - ur 2 3, 
as partial derivatives w.r.t. u and v with the following corners: 

x(0,0)=n, z(l,l) = r 2 , x(l,0)=rg, x(0,l)=r 3 , 



(12) 



(13) 



(14) 



(x(u, v) is our starting surface spanning four straight lines.) For real scalars r and d, we consider two types of 
configurations of the four corners: rulcdi and ruled2. For ruledi we choose 



n = (0,0,0), r 2 = M,0), r s = (0,d,d), n = (r,0,d). 
The mapping from (u, v) to (x, y, z) in this case is 



(15) 



x(u, v) = r (u + v — 2uv) , y{u 1 v)=vd, z(u,v)=ud. (16) 
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For rulecb we choose 

x(0,0)=n, x(l,l)=r 2 , x(0,l) = r 3 , x(l,0) = r 5 . (17) 

The mapping from (it, v) to (x, y, z) in this case is 

x(u,v) = ur 1 y(u,v) = vd, z(u, v) = ud + vd(l — 2m). (18) 

These definitions are such that for r = d the four position vectors ri,r2,r3 and lie at the corners of a regular 
tetrahedron. Fig. [1] and Fig. [2] below are 3D graphs of surfaces called the hyperbolic paraboloids for a choice of 
corners given by (|15[) and (|17|) . . The expression for the mean curvature, calculated using eq. ([3]), of our bilinear 




FIG. 1: The ruledi surface (r — 1, d = 2)with:c,i/ as the horizontal plane and height along z — axis. 
interpolations is 

4r 3 (2u- l)(2v- 1) 
d (d 2 + 2r 2 (2(u - l)u + 2(v - l)v + 1)) 3/2 ' 

for the ruledi and 

4dr(2u-l)(2v-l) 
(d 2 (l - 2v) 2 + 2r 2 (2(u - l)u + 1)) 3/2 ' 



(19) 



(20) 



for the ruled2. 

The mean curvature for our starting surface is zero only for the u = | line and the v = \ line, whereas for a 
minimal surface this should be zero for all values of u and v. Below we describe our effort to improve our surface 
towards being minimal. 



IV. A TECHNIQUE FOR VARIATIONAL IMPROVEMENT 

We want to reduce area of a non-minimal surface x(u, v) using the expectation that reduced value of rms mean 
curvature, denoted by fi, in turn reduces the area A of the surface x(it, v). Our scheme is to reduce the area of a 
surface x(u, v) given by cq. (|10p ~ a special case of eq. ([9|), spanned by a fixed boundary, by obtaining a variationally 
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FIG. 2: The ruled?, surface (r — 1, d = 2) with x, y as the horizontal plane and height along z — axis. 



selected surface X! (u, v) of lesser area. For the variational improvement in surface © , we suggest an ansatz essentially 
consisting of the original surface xq(u,v) of eq. (|10p plus a variational parameter multiplying the numerator of its 
mean curvature. In our notation it becomes 

Xi(u, v, t) = xo(u, v) +tm(u, w)k, (21) 

where t is our variational parameter and 

m(u,v) = uv(l-u)(l-v)H , (22) 

is chosen so that the variation at the boundary curves u = 0,u = l,v = and v = 1 is zero, k is a unit vector 
chosen such that it makes a small angle with the normal to the original surface and Ho, given by ([6]), is numerator 
of the initial mean curvature of the starting surface Xo(it, v). Calling the fundamental magnitudes for xi(u, v) as 
Ei (m, v, t) , F\ [u, v, t) , G\ (it, v, t) , ei (it, v, t) , fi (u, v, t) and g\ (u, v, t), the area A\ of the surface xi (u, v, t) for < 
u < 1 and < v < 1 is given by 

A t = [ [ JExGx -F% dudv. (23) 
Jo Jo v 

We denote the numerator of mean curvature for Xi(u,w) eq. (|21[) as H\(u,v,t). It would have the following familiar 
expression 

H 1 (u,v,t)=E l9l -2F 1 f 1 +G 1 e 1 . (24) 
As Hf(u,v,t) is a polynomial in with real coefficients we rewrite eq. (|24p in the form 

i/ 1 2 ( U , W ,t)=^(/ ll ( W , U ))t\ (25) 

i=0 

Here n turns out to be 10; there being no higher powers of t in the polynomials as it can be seen from the expression 
for E\{u,v,t), Fi(u,v,t) and Gi{u,v,t) which are quadratic in t and e\(u,v,t), f\{u,v,t) and g\{u,v,t) which are 
cubic in t. Integrating (numerically if needed) these coefficients w.r.t. u and v in the range < u,v < 1 we get the 
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following integral for root mean square rms 

Ml (t) = jf Hf(u,v,t) dudv^j = (t* J Y^(hi(u,v)) dudvj . (26) 

The expression in the parentheses on right hand side of above equation turns out to be a polynomial in t of degree n, 
which can be minimized w.r.t. t to find t m i n . The resulting value of t completely specify new surface Xi(m, v). New 
surface Xi(it, v) is expected to have lesser area than that of original surface Xo(it, v). 

In order to see a geometrically meaningful (relative) change in area we calculate the dimension less area by dividing 
the difference of the (original) area of the Coons patch and the variationally decreased area by the original area. 



V. THE TECHNIQUE APPLIED TO HEMIELLIPSOID AND A SURFACE SPANNED BY FOUR 

ARBITRARY LINES 

In this section we apply the technique introduced in the above section[IV]to reduce the area of two types of surfaces. 
In first instance we apply this technique to reduce the area of a non- minimal surface spanning a boundary for which 
the minimal surface is known namely hemiellipsoid eq. (|28p whose boundary is an ellipse lying in a plane and thus 
minimal area in this case is that of the elliptic disc. The reduction in area in this case makes sure the efficiency of 
the algorithm given by cq. (|21l) . In the second example we apply this technique to reduce the area of a bilincarly 
interpolating surface spanned by four boundary lines lying in different planes for which the corresponding minimal 
surface is not known. 



A. HemiEllipsoid-A Surface with Corresponding Known Minimal Surface 



We apply the technique introduced in the section [IV] to the following surface x(u,v) namely hemiellipsoid given by 
cq. (|28p below along with linear blending functions cq. whose boundary is an ellipse. Simpler alternative of the 

above mentioned unit normal N(u,i>), making a small angle with it, in case of hemiellipsoid eq. (|28[) is found to be 

k= (0,0,1). (27) 

A hemiellipsoid 

xo [u, v) — (sin u cos v, b sin u sin v, ccosu). (28) 

with b and c being constants and < u < ^ and < v < 2?t, is a non-minimal surface with its bounding curve an 
ellipse in the ccy-plane; see Fig. [3J In this case we shall treat hemiellipsoid as the initial non-minimal surface and the 
elliptical disc as the minimal surface for the given boundary, namely the ellipse. Thus, eq. (|6]) along with eqs. ([3]) 
and (@| gives 

H = -be sin 3 u (sin 2 u (b 2 cos 2 v + c 2 + sin 2 v) + (b 2 + l) cos 2 u) . (29) 
The root mean square (rms) of beginning curvature given by eq. (|7|) takes the form 

/ pi pi \ !/ 2 

M ° = \J J &2c2sin<5u ( sin2 u ( &2(X)s2 v + ° 2 + sin2 v ) + ( &2 + 1 ) cos2 u ) 2 dudvj . (30) 
The beginning or initial area of the Coons patch given by eq. ((5J takes the form in this case 

A n = J J \J sin 2 u (c 2 sin 2 u (b 2 cos 2 v + sin 2 u) + b 2 cos 2 uj dudv. (31) 

Substituting Hq from eq. (|29|) in cq. (|22jl gives us expression for m(u,v), thus in this case eq. (|21|) becomes 

Xi(u, v) = (sin it cos v, 6 sin it sinw, ccos u — —bct(u — l)u(v — l)v sin 3 u((b 2 — 2c 2 + 1) cos(2u) + 
2(b 2 - 1) sin 2 u cos(2w) + 36 2 + 2c 2 + 3). 
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FIG. 3: A typical hemiellipsoid, initial non-minimal surface of which boundary is an ellipse in the xy — plane 

Finding the fundamental magnitudes E\ (u, v, t) , F\ (it, v, t) , G\ (u, v, t) , e\ (u, v, t) , fx (it, v, t) and gi (u, v, t) for the 
above surface Xi(it,u) eq. (|3"12|) , we can obtain the area A\ using eq. (f2"B")) and H\{u,v,t) using eq. (f2"4"]l and after 
performing the integrations mentioned in eq. (|26[) . the mean square curvature [i\{t) for Xx(w, v) can be calculated. 
These are the similar details as given below for the non-minimal surface spanned by 4— non-coplanar lines. They 
have not been included for this "first instance" but rather included for the "second example" because the formalism 
is well illustrated by this "second example" . Also, that these expressions are too lengthy to be presented. For chosen 
values of b and c we can generate a tabic of their values within the range < u < ^ and < v < 2ir. For our purpose 
we took < b, c < 2 with a step size 0.2 and < u < -| and < v < 2ir , yielding a table of values. Interpolation 
surface for the corresponding minimum values t(b, c) as a function of b and c is given by Fig0] 




FIG. 4: Variation in the parameter t(b, c) as the semi-major and semi-minor axes b and c of the ellipse bounding the hemiellipsoid 
vary. 

In this case the dimensionless decrease p in area for different values of b and c is < p < 23 that may be seen from 
the Fig. [5] 

B. Surface spanned by four arbitrary boundary lines 

Now we apply the technique introduced in the section IIVI to the eq. (fT0|) along with linear blending functions eq. 
for a surface x(it, v) whose boundary is composed of four straight lines connecting four arbitrary corner points 
x(0, 0), x(0, 1), x(l, 0) and x(l,l). For its corners we choose the configuration eq. (fT5|). for a selection of integer 
values of r and d. The results for the configuration (fTF]) have not been included as they are similar to those for the 
configuration (|15j) . We found that the above mentioned simpler alternative of the unit normal N(u,w), making a 



FIG. 5: The dimensionless decrease in area A(b, c) of hemiellipsoid as a function of semi-major and semi-minor axes b and c of 
the ellipse bounding the hemiellipsoid. 



small angle with it, in case of configuration eq. (|T5l) is 

k= (-1,0,0). (33) 
Inserting values of blending functions and boundary points in the eq. (|10[) we find 

Xo(m, v) = (r(u + v — 2uv), vd, ud) , (34) 
with fundamental magnitudes having the expressions as 

E = d 2 + r 2 (l - 2v) 2 , F =r 2 (l-2u)(l-2v), G = d 2 + r 2 (l - 2u) 2 , (35) 

eo = 0, fo = 2d 2 r, and go = 0. (36) 

Thus, eq. ([6]) gives 

H = -4tfV(-l + 2u)(-l + 2v). (37) 
The root mean square (rms) of beginning curvature given by eq. ([7]) takes the form 

^d 2 r 3 

Mo = — ^— ■ (38) 
The beginning or initial area of the Coons patch given by eq. ([5]) takes the form in this case 



A = [ [ d^Jd 2 + 2r 2 (2m 2 ~2u + 2v 2 - 2v + 1) dudv. 
Jo Jo 



(39) 



The scalars r and d can arbitrarily be chosen. Geometrical properties depend only on ratios of lengths, without 
changing the ratio itself and thus without loss of generality d = 1, so that the eq. (|39|) takes the form 



"'O 



A = / / \J Ar 2 u 2 - Ar 2 u + Ar 2 v 2 - 4r 2 u + 2r 2 + 1 dudv. (40) 



(41) 



Substituting H from eq. (|3T|) in eq. (|2"2"j) . we have 

m(u,v) = I6r 3 u 3 v 3 - 24r 3 u 3 v 2 + 8r 3 u 3 v - 24r 3 u 2 v 3 + 36r 3 u 2 v 2 
- \2r 3 u 2 v + 8r 3 uv 3 - 12r 3 uv 2 + 4r 3 uv. 

Using (|4T|) . variationally improved surface eq. ()2ip takes following form 

xi (u, v, t) = [\Qr 3 tu 3 v 3 - 24r 3 tu 3 v 2 + 8r 3 tu 3 v - 24r 3 tu 2 v 3 + 36r 3 tu 2 v 2 - \2r 3 tu 2 v 
+ 8r tuv — 12?* tuv + Ar tuv — 2ruv + ru + rv, v, u). 



(42) 
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(43) 



(44) 



(45) 



Fundamental magnitudes for this variationally improved surface (|42|) are as follows: 

Ei (u, v, t) = (— t(87' 2 (l — u)u(l — v)v(r(l — v) — rv) — 4r(l — u)(l — v)v(r(l — u) — ru)(r(l — v) — rv) 
+ 4ru(l — v)v(r(l — u) — ru)(r(l — v ) — rv)) + r(l — v) — rv) 2 + 1, 

Fi(u, v, t) = (— i(8r 2 (l — u)u(l — v)v(r(l — u) — ru) — 4r(l — u)u(l — w)(r(l — u) — ru) 
(r(l — v) — to) + 4r(l — u)uv(r(l — u) — ru)(r(l — v) — rv)) + r(l — u) 
— ru)(—t(8r 2 (l — u)u(l — v)v(r(l — v ) — rv) — 4r(l — u)(l — v)v(r(l — u) — ru) 
(r(l — v) — rv) + 4ru(l — v)v(r{l — u) — ru)(r(l — v) — rv)) + r(l — v) — rv), 

Gi(u,v,t) = (-t(8r 2 (l - u)u(l - v)v(r(l - u) - ru) - Ar{\ - u)u{\ - v)(r(l - u) - ru) 

(r(l — v) — rv) + 4r(l — u)uv(r(l — u) — ru)(r(l — v) — rv)) + r(l — u) — ru) 2 + 1, 

ei(u, v, t) = i(-96rW + 144r 3 wu 2 - A8r 3 uv + 48rV - 72rV + 24r 3 v), (46) 

fl(u,v,t) = t(-U4r 3 u 2 v 2 + lUr 3 u 2 v - 24r 3 w 2 + 144r 3 uu 2 - 144r 3 wv + 24r 3 u - 24rV + 24r 3 u - 4r 3 ) + 2r, (47) 
and 

gi (u, v, t) = t(-96r 3 u 3 v + 48r 3 u 3 + 144r 3 u 2 w - 72r 3 it 2 - 48r 3 uu + 24r 3 u). (48) 

Inserting these values of fundamental magnitudes in eq. (|24| we find the expression for Hi(u, v, t) of surface (j42]l as 

Hi(u,v,t) = [-4r 3 (2w- l)(2v - 1)] + [8r 3 (2u - l){2v - l){r 2 (u 2 (6v - 5){6v - 1) + u(-36(v - l)v - 5) + 5(v - 1) 
v + 1) - 3(u 2 + v 2 ) + 3{u + v))]t + [-32r 7 (2u - l)(2v - l)(6u 4 (2(w - l)v(18(v - l)v + 5) + 1) - 12u 3 
(2(v - l)v{18{v - l)v + 5) + 1) + u 2 (2{v - l)v(138(v - l)v + 37) + 7) + u(-2(v - l)v(30(v - 1) 
v + 7) - 1) + (v - l)v(6(v - l)v + l))]t 2 + [u(2u - l)(v - l)v(2v - l)(12u 4 (3(w - l)v(l2(v - 1) 
v + 5) + 2) - 24w 3 (3(w - l)v(12(v - l)v + 5) + 2) + 3u 2 (12(t; - l)v{17{v - l)v + 7) + 11) 
- 9u(A(v - l)v(5{v - l)v + 2) + 1) + 3(v - l)v{8{v - l)v + 3) + 1)]< 3 . 

(49) 

After performing the integrations mentioned in eq. (|26|). the mean square curvature ^ 2 (t) for Xi(u,w) becomes 



^2048r 18 \ 6 /190464r 16 \ 5 | ^512r 12 (l53r 2 + 77) t ( 256r 10 (7r 2 + 3) 
7275/ t + V 25050025 ) * 

'32r 6 (29r 4 + 98r 2 + 119)' 



^ 2277275 ) 1 + ^ 25050025 ) 1 + [ 444675 j* + [ 3675 J* 

7 (50) 



1225 



which may be minimized for t for every fixed value of r. Fig. [6] represents this minimizing value of t m in as the 
numerical function of r. 

We find the variationally improved surface Xi (u, v) eq. (|2ip and its area as given by eq. (|23[) for each t m i n for the 
corresponding r. For a selection of r values for < r < 2 with step size 0.001, the dimension less decrease in area of 
surface x (u, v) of eq. (fTU)) can be seen in the Fig. [7] and interpolating curve of the same is provided in Fig. [5J 

VI. CONCLUSIONS 

We have discussed a technique to reduce the area of a surface x(it, v) eq. ([9]) obtaining variationally improved 
surface Xi(u,w) of eq. ()21j) . This algorithm is first applied to a non- minimal surface spanning a boundary for which 
minimal surface is known namely hemicllipsoid eq. (1281) . The dimensionless decrease p in the- area of the hemicllipsoid 
eq. (|28[) for different values of b and c is < p < 23 (see Fig. [5]) depending upon how much it is far from the minimal 
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FIG. 6: The variation in parameter t(r) depends on the variation of real scalar r, for the skew quadrilaterals ruledi bounded 
by four arbitrary straight lines connecting four corners x(0,0),x(0, l),a;(l,0) and x(l, 1); 
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FIG. 7: The dimensionless decrease in area A(r) as a function of r enclosed by a skew quadrilateral ruledi bounded by four 
arbitrary straight lines connecting four corners x(0,0),x(0, l),x(l,0) and x(l, 1). 
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FIG. 8: The dimensionless decrease in area A(r) as a function of r enclosed by a skew quadrilateral ruledi bounded by four 
arbitrary straight lines connecting four corners x(0,0),x(0, l),x(l,0) and x(l, 1). 

surface, namely the elliptic disk. This shows our algorithm cq. (|21[) can significantly reduce area of surface that is 
far from being minimal. After noting this effectiveness, we applied this technique to reduce the area of a surface of 
x(u, v) eq. (fTUj) bilinearly spanned by four non-planar boundary lines, a special case of Coons patch eq. ©, along 
with the configuration eq. (fT5)) , for a selection of r values for < r < 2 with step size 0.001. This gave us a much 
lesser ( in the range to 0.80) dimensionless decrease in less area of surface x (u, v) of eq. (fT0|) , as seen in the Fig. [7] 
or Fig. |S1 This suggests that ruledi is already a near minimal surface. 
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